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(ABSTRACT) 


Studies on two adhesives (Metibond 1113 and' 1113-2) identified 
as having applications in the bonding of composite materials are pre- 
sented. Constitutive equations capable of describing changes in 
material behavior with strain rate are derived from various theoretical' 
approaches. It is shown that certain unique relationships exist be- 
tween these approaches. It is also shown that the -constitutive equa-- 
tion derived from mechanical models can be used for creep and relaxation 
loading. 

Modifications to the constitutive equations for constant strain 
rate loading are proposed. Nonlinear constitutive equations are 
derived-using a nonlinear perturbation technique in conjunction with a 
modified Bingham model. Using the modified Bingham model, constitutive 
equations describing loading-unloading behavior are also proposed. 

The stress-strain behavior of the adhesives is shown to be 
significantly rate dependent.' Further it is shown that a rate dependent 
stress-whitening (crazing) phenomenon occurs prior to the maximum 
stress. A region of elastic behavior, a rate and time dependent region, 
and a region of perfectly plastic flow are identified in the stress- 
strain behavior. Information regarding variations of Poisson's ratio 
with rate and time is also presented. 

The elastic limit stress and strain, and maximum stress are shown 
to be rate dependent and agree well with an empirical equation pro- 


i 



posed by Ludwik. Analytical predictions based on modified Ramberg- 
Osgood equations are shown to' agree well with experimental s'tress- 
strain-strain rate results. It is shown that the coefficients of 
these equations are different before and after stress-whitening due 
changes in the properties of the adhesives. Analytical predictions 
based on the modified Bingham model are shown to agree well with the 
constant strain rate results. It is also shown that the nonlinear 
model indicates that the coefficients of the modified Bingham model 
may vary due to the change in material properties before and after 
s t res s -wh i ten j n g . 

A creep to failure phenomenon is shown to exist and is correlated 
with a delayed yield equation proposed by Crochet. Loading-unloading 
results are presented and are shown to correlate well v/ith the proposed 
form of the loading-unloading equations for the modified Bingham model. 

Experimental results obtained for relaxation tests above and 
below the glass transition temperature are presented. It is shown that 
the adhesives obey the time-temperature superposition principle. 
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Chapter 1 
INTRODUCTION 

Adhesives are currently being used in a wide variety of structur- . 
al applications. This is especially true for structures made using 
advanced composites. Adhesive bonding has proven to be a viable method 
of assembling composite joints, however the overall performance of an 
adhesively bonded structural composite cannot be predicted accurately 
unless the role of the adhesive joint is fully understood. 

It has been pointed out that the process of adhesion can be 
divided into two types; mechanical adhesion and specific adhesion [1].* 
The former describes the process of the adhesive solidifying in the 
pores of the two adherent surfaces and the latter describes the 
process of attractive intermolecular forces between molecules of even 
the smoothest solids. For mechanical-adhesion. Bowers and Zisman [1] 
state that if ideal conditions (i.e., complete wetting and freedom from 
the formation of gas pockets arid inclusions) prevail, the joint must 
fail in cohesion rather than adhesion. That is, failure is in the bulk 
phase rather than at the adherend/adhesive Interface. Since bulk or 
cohesive failure is a likely failure mode, the present study is an 
effort to identify bulk or cohesive stress-strain, strain-rate, creep, 
relaxation, yield and/or failure properties of typical structural 

*Numbers in brackets [ ] refer to the references listed in the 
Bibliography. 
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adhesives. This is in contrast to the work of others regarding the 
viscoelastic behavior of adhesives [2] where the adhesives were tested 
in the bonded state. 

A complete analysis of the stress distribution in adhesively 
bonded structural composite joints must include rate and time dependent 
material properties when these are significant. Adhesive materials, 
which are frequently molecular high polymers, generally exhibit signifi- 
cant rate and time dependence that must be considered in a reliable 
failure analysis. The phenomenological behavior of the adhesives in 
bulk form is considered to be of prime importance in this study. The 
extent to which bulk properties can be related to properties in the 
bonded state is as yet undetermined. However, it is reasonable to 
assume that if adhesive materials exhibit rate and time dependent 
properties in bulk form, time effects will also be important in the 
bonded state. In addition, this study may give insight into rate and 
time dependent behavior of some matrix materials used in advaqced 
composites as they are frequently similar in composition to the adhe- 
sives investigated herein. 

With the advent of computer oriented analytical techniques to 
study the behavior of complex structural problems, it is now possible 
to observe the effects of the plastic and viscoelastic properties of 
the constituent materials. The performance of a structural system is 
predicted from the mechanical properties of the materials employed. 

If the environment causes strain rates above the static rate, the re- 
sulting effects on modulus, strength, and ductility must be understood. 



An analytical description of mechanical behavior by a suitable 
constitutive equation accounting for property changes with strain, 
rate becomes necessary for accurate response calculations. 

The materials under investigation .in this study are Narmco 
Whittaker's Metlbond lllS and 1113-2 adhesives. These adhesives are 
currently being, used in composite- bonding applications. Initial 
investigations indicated that the adhesives exhibited strain rate 
effects with respect to elastic limit stress and maximum stress. For 
a given strain rate, a region of linear elas.ticity and a region of 
inelastic behavior followed by perfectly plastic flow at the maximum 
stress to failure was observed for both adhesives (these results can 
be found in Reference 3). It was also observed that a "stress-whiteni 
phenomenon occurred prior to failure. The constitutive equations 
describing these materials should include as many of these rate and 
time dependent properties as possible. 

The second chapter of this dissertation is devoted to selecting 
the final constitutive forms to be further analyzed. A review of 
constitutive equations derived from various theoretical approaches is 
presented. These are deformation theories,' incremental theories, and 
mechanical models. It is shown that certain unique relationships 
exist between these approaches. The final section of this chapter 
concerns the use of the constitutive equation derived from mechanical 
models for creep and relaxation loading. 

Chapter three outlines proposed modifications to the constitutive 
equations. A nonlinear perturbation technique developed by Davis [4] 
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used in conjunction with the constitutive equation derived from 
mechanical models is reviewed. The use of this same mechanical model 
for loading-unloading considerations is also reviewed. 

Chapter four describes the experimental considerations. The 
materials selected for study and specimen features are reviewed. A 
description is also given of the experimental apparatus used for the 
various types of tests conducted. 

The fifth chapter is concerned with' the constant strain rate 
results and the ability of the constitutive equations selected to 
model this behavior. It will be shown that during the constant strain 
rate testing a "stress-whitening" or "crazing" phenomenon occurs prior 
to the maximum stress.. The effect of this phenomenon on the 
constitutive equations is reviewed. It is also shown that the 
perturbation technique developed' by Davis yields significant. information 
relative to the behavior of the adhesives above and below the 
"stress-whitening" point.' 

Chapter six concerns the ability of the constitutive equation 
derived from mechanical models to predict the behavior of Metlbond 
1113 and 1113-2 for creep and relaxation as well as other more complex 
loading histories. This chapter includes comparisons between 
theoretical predictions and experimental results for two different 
loading histories. 

Since bonded .composite structures are often subjected to high 
temperature environments, the behavior of the adhesive at high 
temperatures is also important. Chapter seven presents the experimental 
results obtained for relaxation tests above and below the glass 
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transition temperature. Theoretical considerations relative to the 
time-temperature superposition principle are presented. 

The final chapter presents the conclusions and indicates areas 
for future research. 



Chapter. 2 

CONSTITUTIVE EQUATIONS 

This chapter reviews the development of constitutive equations 
which have been used previously to model rate and time dependent 
material behavior. These equations were derived from several ap- 
proaches; deformation theories, incremental theories, and mechanical 
models. 

It is shown that certain unique relationships exist between the 
final form of the constitutive equations derived from incremental 
theories and that derived from mechanical models. It is also shown 
that constitutive equations proposed on the basis of experimental ob- 
servations by early investigators can be derived from recent theories 
and mechanical models. The final section of this chapter concerns the 
constitutive equation derived from mechanical models for creep and 
relaxation loading. 

2.1 Terminology 

During the development of the constitutive equations certain 
terms such as viscoelastic, viscoplastic, yield stress, and work- 
hardening are used. Therefore, before proceeding into the various 
approaches these terms are discussed. 

The stress-strain diagram shown in Fig. 1(a) is representative 
of the materials studied herein. It is assumed that this response 


6 



7 


would be obtained if the material were loaded- or unloaded at a constant 
rate. Several different types of analytical approaches to the charac- 
terization of such stress strain response are given. ' In each, the 
parameters e and Y have different meaning and the general approach to 
the representation of the region beyond e is different in each case. 

. The modified Bingham model or linear, viscoelastic-plastic mechan- 
ical model of F-ig. 1(b) is used in which e is an elastic limit stress 
and Y is a yield stress. The material is linear elastic below 9, 
linearly viscoelastic between 0 and Y and perfectly plastic above Y. 

A non-linear viscoelastic modified Bingham model is. also dis- 
cussed. For this model the material is elastic below 6, nonlinearly 
viscoelastic between 0 and Y and perfectly plastic above Y. 

In all other instances (i.e., the Ramberg-Osgood approach and 
much of the discussion of literature), the material is assumed to be 
elastic below e, strain hardening plastic (including time effects) 
between e arid Y and perfectly plastic beyond Y. 

Without exception all comparisons between theory and experiment 
are made using either the linear viscoelastic-plastic modified Bingham 
model or the non-linear modified Ramberg-Osgood model. 

The term viscoplastic will refer to rate effects occuring in the. 
plastic region. Perzyna [5] defined an "elastic-viscoplastic" material 

t 1 

as having rate effects in both the elastic and plastic regions. In 
contrast an "elastic/viscoplastic" material' was defined as having rate 
effects in the plastic region only. This terminology has been adopted 
for this study. 
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Figure 1. Stress-Strain Curve and .Modi fi ed Bingham Model 
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2,2 Defomation Theories 

In 1924, Hencky [6] proposed stress-stress relations whereby the 
plastic strains are functions of the current state of stress and are 
independent of the loading history. Such theories are called total or 
deformation theories. Hencky's relations, as given by Kachanov' [7] 
are, 

j " 9f“ ^i j (2.2.1) 


where K = E/3(l - 2v) is the bulk modulus and E and v are the elastic 

modulus and Poisson's ratio, respectively. The term T is a scalar 

function of the invariants of the stress tensor as indicated by Hill 

[8]. The parameters a-jj, and represent the stress, strain, 

and deviatoric stress tensors, respectively. 

For f = constant = 1/2G, where_G_= E/2(l +'v) is the shear 

modulus, Eq. (2.2.1) reduces to the elastic stress-strain relations. 

Ho\yever, If f = 1/2G + iji, the strain components are then the sum of 

E P 

elastic components, e-jj , and plastic components, , and Eq. 

(2.2,1) becomes. 


^mm ^ 
^ij " 9K 


S - - 

2G 


+ S 


ij 


( 2 . 2 . 2 ) 


or • 


where 


®ij ^1j ®1j 

E „ 5m j j. ^ 
^ij 9K ^ij 2G 


(2.2.3) 


(2.2.4) 
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^ S^-j (2,2.5) 

and is a scalar function of the invariants of the stress tensor. 

For the case of uniaxial tension, Eq. (2.2.2) reduces to, 

e = ^+<j)a (2.2.6) 

2 

where the subscripts have been dropped and 4) = :r 4). 

The function 4> in Eq. (2.2.6) is often determined empirically in 
order to fit experimental results. 

2.2.1 Ramberg-Osgood Equation 

From uniaxial tension and compression results on various metal 
alloys, Ramberg and Osgood [9] in 1943 proposed an analytical expres- 
sion of the form. 


e = I + Ka" (2.2.7) 

where K and n are material constants. It is evident that Eq. (2.2.7) 
may be considered as a special case of Eq. (2.2.6) since 4> is a 
function of 0 . 

2.2.2 Modified Ramberg-Osgood Equation 

In 1960 McLellan [10] modified Eq. (2.2.7) for rate dependent 
materials by observing that for many materials n is invariant with 
strain rate and K and E are simple functions of the strain rate such 
that, 


cf'-i 'd. 
E(e) = ce 


( 2 . 2 . 8 ) 



n 


and 

K(e) = • (2.2.9) 

where a, b, c, and d are material constants. Variations in E, K, and 
n (if such exi.sts for n) indicate different types of material sensi- 
tivity to strain rate. Changes in E denote viscoelastic behavior. 

Both parameters K and n of the plastic strain term represent work- 
hardening characteristics and variations in these parameters with 
strain rate can be considered viscoplastic phenomenon. Changes in K 
show that the relative magnitude of strength is affected by strain 
rate, whereas changes in the plastic flow behavior are indicated by 
variations in n. 

Substituting Eqs. (2.2.‘8) and (2.2.9) .into Eq. (2.2.7) McLellan 
obtained the form, 

e.= -?a-+aA" (2.2.10) 

McLellan and Eichenberger [11] in 1967 used an equation of this form 
to predict the behavior of aluminum in uniaxial compression at various 
strain rates. McLellan [12], also in 1967, was able to modify Eq. 
(2.2.10) in order to predict the behavior of many metals, various 
phenolics, and plexiglass at various strain rates and temperatures. 

In 1971 Zabora, et al. [13] used a form of Eq. (2.2.10) to describe the 
shear behavior of structural adhesives in the bonded state. 

Eq. (2.2.7) and its respective rate form, Eq. (2.2.10), were 
selected for further consideration in this study. 


■ i 
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2.3 Incremental Theories 

In the incremental theory of plasticity the strains are in 
general not uniquely determined by the stresses, but depend on the 
entire history of loading. The equations describing plastic flow can- 
not in general be finite relations between the components of stress 

and strain, but must be incremental relations. 

\ 

The equations of incremental theory establish relations between 
stress and infinitesimal increments of strain. Equations of this form 
were proposed in 1925 by Prandtl [14] and later in 1930 by Reuss [15] 
and are known as the Prandt-Reuss equations (textbook presentation can 
be found in Mendel son [16]). 

Reuss assumed that the plastic strain increments are at any 
instant of loading above yield proportional to the instantaneous stress 
deviation, 

de^j = dXS-jj (2.3.1) 

where dA is a non-negative parameter which may vary throughout the 

loading history. To determine the actual magnitudes of the plastic 

strain increments a yield criterion is required. 

In incremental theory the total increments of deviatoric strain, 

• E 

Ae^-j', are the sum of elastic components, Ae^*j , and plastic components, 

p 

Ae^-j . In rate form (i.e., dividing by At and taking the limit as At 
approaches zero) this becomes. 


• E P 
®ij ®ij ®ij 


(2.3.2) 
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. E 


'ij 


For an elastic/viscoplastic material, the elastic strain rate, 
, is given by. 



(2.3.3) 


p 

The plastic strain rate, e. . represents combined viscous and plastic 
effects and thus is often called the "viscoplastic component." Possi- 
ble constitutive forms for the viscoplastic component were given by 
Hohenemser and Prager [17] in 1932 in the form. 


where 



2k<F> 




(2.3.4) 


(2.3.5) 


is the yield function, k is the yield stress in simple shear, y is the 
viscosity coefficient, J2 represents the second invariant of the 
stress deviation given by, 


*1 2 2 2 
•^2 " “ ^ 22 ) ^*^22 ■ ^ 33 ) (^^33 ” '^ ll ) 


+ 6(a-|2^ + + <^31^)] 


( 2 . 3 . 6 ) 


and the symbol <F> is defined, 


0 for F _< 0 

F for F > 0 


(2.3.7) 


The yield function, Eq. (2.3.5), implies von Mises yield 

1 /2 

criteria. That is, yield is said to occur when Jg = k. The total 
strain rate below yield is the elastic strain rate. Above yield, the 
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total strain rate is composed of elastic and viscoplastic components. 
Combining Eqs. (2.3.2) through (2.3.5) the following relations are 
obtained, 


; - hi. 

®ij ■ 26 

S-. l-kjl/2 

®i j, = 2 ^ + 2^^ ^i j 


1 V2 . 
k 

1 u 

^ k 


(2.3.8) 


Equations of this form were studied by Freudenthal [18] in 1958. 

In 1963 Perzyna [19] generalized Eqs. (2.3.8) by replacing 
2k<F> by y°^(F), where y° denotes a material' constant, F is the yield 
function (not necessarily given by Eq. (2.3.5)), and the function $ 
satisfies' the conditions. 


$(F.) = 0 for F _< 0 

$(F) M for F > 0 


Perzyna 's generalized form was. 


, = ^ 

'ij 26 


F < 0 


®ij = I" > 0 


3F 


(2.3.9) 


(2.3.10) 


where y = y°/2u is a material constant. 

The relations in Eqs. (2.3.10) involve the assumption that the 
viscoplastic component be a function of the over-stress above the 
yield condition. This assumption is the same as that introduced by 
Malvern [20] earlier in 1951 from experimental investigations of one- 


i 
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dimensional problems. 

2.3.1 Sokolovsky's. Equation 

During the same era as Malvern's work, Sokolovsky [21] (1948) 
determined from experimental observations that the viscoplastic com- 
ponent can be written as a function of the over-stress above the 
elastic limit, 9(0 - e). Sokolovsky's equations had the form, 

• 0 

e = ^ . 0 0 

(2.3.11) 

o 1 

£ 9(c^“0) 0^0 


It is interesting to- note that the form of Sokolovsky's equations 
can be determined from that proposed by' Perzyn a, Eqs. (2.3.10), if the 
yield condition is taken as Eq. (2.3.5). The function $(F) is then 
written as. 


fo 


• Hn = 


IF 


for F _< 0_ 
for F > 0 


or 1 k 

or > k 


(2.3.12) 


and Eqs. (2.3.10) become. 


s. • 

; ''J 

®ij ■ 2G 


®i j = # + T 


1/2 


- 1 


^ij 




1 V2 . 
Jo < k 


1 1/2 . 
J2 ^ k 


(2.3.13) 


For one-dimensional loading (and dropping the subscripts), Eqs. 
(2.3.13) become. 
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E 



cr < /3k. 


£. = 




(a - /3k) 


a > /Sk 


(2.3.14) 


For yon Mises criteria, the yield stress in pure shear is 1//3 times 
the yield stress, in simple tension. That is, /3k = e, and hence Eqs. 

(2.3.14) have the form. 


e = 


.e 



o <^'e 
cf > 6 


(2.3.15) 


— 2 

where y = w v^y is a material constant. These are essentially 

( 

Sokolovsky's equations. 

The form of Eqs. (2.3.15) was selected for further consideration 
in, this study. The final equations investigated are actually derived 
in the following section using a mechanical model. However, Eqs. 

(2.3,15) yielded information relative to the coefficient of the over- 
stress term which was not readily recognizable in the mechanical model. 
The elastic limit or yield stress of the adhesives was found to be 
rate dependent (i.e., e = e(e)). This then indicates that the cor 
efficient of the over-stress term in Eqs. (2.3.15) may be rate 
dependent. 


2.4 Mechanical Models 

An alternate approach to the development of suitable constitutive 
forms is the use of discrete mechanical models. Recent investigators 
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[ 22 , 23 , 24] have attempted to use such models to describe the. rate 
dependency observed in various polymers. 


2.4.1 Modified Bingham Model 

In 1972, Brinson [23] proposed a modified Bingham model. Fig. 
1(b), to describe the rate behavior of various polymeric materials. 
The constitutive equations for this model are. 



a < 0 


e 



e < a < Y 


(2.4.1) 


Once the maximum stress, Y,- is reached, perfectly plastic flow to 
failure is assumed. It can be readily observed that Eqs. (2.4.1) are 
similar to those proposed by Sokolovsky, Eqs. (2.3.11), and that 
derived from the work of Perzyna, Eqs. (2.3.15).- 

Using Eqs. (2.4.1) to model the behavior of polycarbonate at 
different strain rates, Brinson observed that the viscosity coefficient 
had to vary with strain rate in order to accurately represent the 
experimental results. Brinson and DasGupta [24] in 1973 reported the 
same phenomenon using a series of Bingham elements. 

It could not be ascertained directly .from Eqs. (2;4.1) that y 
(the coefficient of the over-stress term) was rate dependent. However, 
the analogous form of these equations, Eqs. (2.3.15), indicated that 
this coefficient may be rate dependent. Thus, insight into rate 
dependency not readily seen using mechanical models can be observed 
from incremental theory as derived for an elastic/viscoplastic material 


'I I ‘ ’ 11 H 


M ' 
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This rate dependency in p is later shown (Chapter 5) to be an inherent 
property of rate-sensitive plastic materials. 


2.4. 1.1 Constant Strain Rate Solution 

For constant strain rate (e = R = constant), Eqs. (2.4.1) are 
solved according to the initial condition, 

■ “ ^0^ ~ ®o (2.4.2) 

where and tg are strain and time, respectively, at which the 
elastic limit is reached. The solution is. 


0(e) = eE 


0 < 9 


0 (e) = 0 + pR{l - e"^?^ “ e < o < Y 


(2.4.3) 


2. 4. 1.2 Creep and Relaxation Solutions 

The constitutive equations reviewed have been concerned with 
modeling the effects of strain rate on material behavior. It would be 
ideal if these equations modeled material behavior for other loading 
conditions. 

Equations derived from deformation theory would not be suited to - 
this task since in deformation theory the plastic strains are a. 
function of the current state' of stress and are independent of the • 
loading history. Equations derived from the incremental and mechanical 
model approaches are written in incremental form and thus can take into 
account the loading history. Therefore, the modified Bingham, model 
was considered further for other than monoton i call y increasing, loads. 
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For creep loading using the modified Bingham model. 


0(t) = 0QH(t) > 0 


(2.4.4) 


where H(t) represents the unit step function and Oq is the level of 
constant stress. Eqs. (2.4.1) then become. 


e = 


00 - e 


0 > 0 


(2.4.5) 


The initial condition is given as, 


:(t = 0) = ^ 


(2.4.6) 


That is, the initial response is assumed to be instantaneous and 
elastic, and the solution to Eq. (2.4.6) is written as, 


:(t) = 


oo - e 




(2.4.7) 


For relaxation loading, 

e(t) = ioH(t) > Sq (2.4.8) 

where Fq is the level of constant strain and Eq elastic limit 

strain. Eqs. (2.4.1) then becomes 

|.+ i (a - 0 ) = 0 0 > e (2.4.9) 

■ E p ' 

subject to the initial condition, 

o(t = 0) = EqE (2.4.10) 

That is, the initial response is assumed to be instantaneous and 
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elastic. 


The solution to Eq. (2.4.9) is written as, 

, , , -Et/y 

a(t) = e + (sqE - 9 )e 


(2.4.11) 



Chapter 3 

PROPOSED CONSTITUTIVE EQUATIONS 


3.1 Modified Ratnberg-Osgood Models 

The plastic strain term, Kcr*^, in Eq, (2.2,7) is a function of 
the current state of stress. It is proposed that the plastic strain 
be a function of the over-stress (a - e). The proposed form is. 


^ ~ E . . 


cr < e 


e = ^, + K(ct “-6)*^ ■ ' a •> 6 


(3.T.1) 


where K arid -n are material constants. 

Following a procedure similar to that -of McLellan [10]', section 
2.1.2, Eqs. (3.1.T) were modified for rate effects. The final form is. 


e = -7d ^ ^ 

■ (3.1.2) 

e = “TT t ae^(a - 9)*^ cr > 0 
ce*-* 


where a, b, c, and d are material constants. Variations in K and n 
with strain rate again denote work-hardening characteristics. 


3.2 Nonlinear Model 

The following section is a review of a nonlinear perturbation 
technique used in conjunction with the modified Bingham model in an 
attempt to describe the rate effect observed in the viscosity coef- 
ficient. 
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111! 1 U’lt If. (I" 1 
n;‘‘t' ' 1.. !' ■ 

The basic hypothesis in mechanical models of linear visco- 
elasticity is that the elastic and viscous coefficients — parameters 
representing the spring and dashpot, respectively— are independent of 
the strain and strain rate. In this section this requirement is 
relaxed by assuming a small term depending on the elastic strain is 
added to the elastic coefficient, and. a small term depending on the 
strain rate in the viscous element is added to the viscous coefficient. 

Using an hypothesis proposed by Davis [4] in 1964 for'first 
order nonlinearity, the nonlinear modulus and viscosity coefficient 
can be written as, 


^ ^o/ teEoEe , 

Eo . 1/ 

+ Yu £ij 


(3.2.1) 


where and Vq represent the linear modulus and viscosity coefficient 
of the spring and dashpot, respectively, and X = Eg/yo represents the 
inverse of the. linear relaxation time, t. The parameters and e 
are the strains in the spring and dashpot, respectively. The coef- 
ficients and Y^ are dimensionless constants which are a measure of 
the nonlinearity. The conditions of first-order nonlinearity are. 


Iy^e^I « 1 
t|y e I « 1 


(3.2.2) 


1/ 

For elasttc/viscoplastic materials the term L could be 
represented by the viscoplastic component However,- the above termi- 
nology' is used in this- development in order-to distinguish between the 
elastic and viscous elements. 
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and higher order terms in ye^e and ane neglected. 

The nonlinear constitutive equations for the modified Bingham 
model using the relationships in Eqs. (3.2.1) are. 


e 




ao 



1_ 

I'o 


(a - e) 



(a - e)^ 

VoEo 


CT £ 9 

(3.2.3) 

0 < a < Y 


Eqs. (3.2.3) are derived in Appendix A. It can also be seen that when 
Te = = 0 the linear form Eqs. (2.4.1) is obtained. 

The perturbation technique proposed by Davis can now be used to 
solve Eqs. (3.2.3) for the condition of constant strain rate. Expand- 
ing a in a power series in terms of and Yy of the form. 


00 p CO 

a = Oq ^ ^ TySp (3.2.4) 

n=l n=l ^ 

and limiting this expression to first-order nonlinearities gives, 

a = Oq + Tpa-] + Yy$i (3.2.5) 

Eq. (3.2.5) is substituted into Eqs. (3.2.3) and the expressions 
for aQ, a-j , and S-| determined (see Appendix A). Postulating, that non- 
linear effects are present in the dashpot only (i.e., Yp = 0) the 
following nonlinear constitutive equation is realized for the modified 
Bingham model , 



24 


a(e) = Eq£ g 5 9 

= 8 + PoR(l - e ^o)/^'oR) + y^(2EqRt(£ - Eq) 


(3.2.6) 


+ EoR^t^O 


^-2E(e - eo)/Vo'^) 


9 < a < Y 


Eqs. (3.2.6) are also derived in Appendix A. 

As stated earlier, the linear coefficient pq was found to be 
dependent on the total strain rate e.- The nonlinear hypothesis, Eqs. 
(3.2.1), indicates that this linear value may vary with respect to 
the strain rate in the viscous element. Thus,' the nonlinear model 
does not mathematically indicate that pq may be dependent on e, but 
rather that it may change with respect to However, it is shown 
(Chapter 5) that the nonlinear solution indicates a change in the 
behavior of the adhesives above and below the stress-whitening point. 


3.3 Loading-Unloading Considerations 

The work-hardening characteristics of the adhesives were found 
to be rate dependent. That is, the maximum stress and plastic flow 
behavior varied with strain rate. However, the elastic limit stress, 
the point of initiation of viscous and plastic effects, remained 
constant during loading and unloading cycles at a given strain rate. 
The adhesives were also observed to unload elastically. 

For the considerations above, the function 'i'(F) in Eq. (2.3.9) 


can be modified as, 



r 




0 


^(F) = 


F 


for F £ 0 
or F < 0 
for F > 0 
and F > 0 


(3.3.1) 


For one-dimensional loading and von Mises yield criteria, these con- 
ditions become. 


0 
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for a £ 0 
or a < 0 
for a > 9 
and 0 > 0 


(3.3.2) 


The conditions given by -Eq. 
modified Bingham model as, 



(3.3.2) can be incorporated into the 

for a 5 . 0 
or 0 < 0 

(3.3.3) 

for 0 < 0 < Y 
and 0 > 0 


Eqs. (3.3.3) constitute the proposed form of the modified Bingham 
model for loading-unloading considerations-. 
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Chapter 4 

EXPERIMENTAL CONSIDERATIONS 

4.1 Materials and Specimen Features 

The' two adhesives Investigated in this study were Narmco Whit- 
taker’s Metlbond 11T3 and 1113-2. Metlbond 1113 is a 100% solids, 
modified epoxy film with a synthetic carrier cloth. Metlbond 1113-2 
is the identical film without the carrier cloth. Tests were run on 
both materials to ascertain the behavior of the modified epoxy in both 
supported and unsupported form, that is, to establish the effect of 
the carrier cloth by comparison. 

All tests were conducted on bulk specimens such as that shown in 
Fig. 2, Load was introduced through steel plates bolted to the 
specimen and aligned using a special alignment fixture. The specimens 
were cut from -v 0.140 in. ('u 0.355 cm.) thick sheets which were laid 
up using fourteen plies of film and cured according to manufacturer's 
specifications [28] in a platen press at NASA-Langley Research Center. 
These sheets were 12 in. (a. 30.5 cm.) square. The width of each 
specimen was -\i 0.5 in. {0/1.27 cm.) and the gauge length was o/ 3.25 in, 
(o. 8.26 cm.). Initially six sheets of each adhesive were fabricated 
and either six or seven specimens were cut from each sheet. These 
specimens were used to conduct constant strain rate, creep (constant 
stress), and relaxation (constant strain) tests. Six more sheets of 
each adhesive were then fabricated from a second batch of material. 

The majority of specimens from these sheets were used for the loading- 
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Figure 2. Specimen Mounted for Testing 
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unloading tests and for the high temperature testing. Over 90 speci- 
mens were tested in all. Controlled environmental conditions of 
approximately 72° F and 75% relative humidity were used in all tests. 

4.2 Experimental Procedure 
4.2.1 Constant Strain Rate Tests 

Constant strain rate tests were performed on an Instron testing 
machine using constant head rates ranging from 0.002 in. /min. (0.00508 
cm. /min.) to 2 in. /min. (5.08 cm./min.). The experimental apparatus 
for these experiments is shown in Fig. 3. Each specimen was instru- 
mented with an Instron model G-51-13 extensometer (Fig. 2) and both 
longitudinal and transverse electrical resistance foil type strain 
gauges (M.M. EP-08-125BB-120 LE). The extensometer was used to ascer- 
tain if heating or reinforcement of the electrical gauges affected 
strain gauge results. There was good agreement between both methods 
for nearly the full range of strain. A small difference in strain 
values was observed after local yielding (i.e., the formation of craz- 
ing or microcracks), with the electrical gauges giving slightly lower 
results. This behavior has been attributed to the slight reinforcing 
effect of the strain gauge. The signals from the extensometer and the 
foil gauges were amplified using Vishay (Model BAM-1) bridge amplifiers 
and were recorded on a Hewlett Packard (7100 B) dual channel strip 
chart recorder. The bridge amplifiers were operated at a reduced 
voltage which allowed recording of strains up to 10% (0.10 m/m). 

The reduced voltage also lowered strain gauge heating effects. 
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Figure 3. Constant Strain Rate Testing Apparatus 
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To establish the repeatability of results, two specimens from 
one sheet and one specimen from a second sheet of each adhesive were 
tested at the same strain rate. In all cases, differences of less 
than 2% were found in the stress-strain behavior below failure. The 
variation in the failure stresses varied similarly, but failure 
strains were more random. As a general rule, results were obtained on 
the bases of a single test. However, a second specimen was tested 
whenever a specimen failed prematurely due to defects in the specimen. 
Such defects were occasionally present in the form of internal bubbles 
which formed during the curing process. 

4.2.2 Creep Tests 

Creep tests were performed using a pneumatic testing machine 
(manufactured by Allied Research Associates) which was capable of 
loading at a rate of approximately 20 in. /min. (50.8 cm. /min.) and 
then maintaining a constant load. This apparatus is shown in Fig. 4. 
The load level was monitored using a Baldwin SR-4 (Type U-1 ) load cell 
with a Baldwin (Type L) strain indicator. The load was constant 
throughout the duration of each test. Since the reduction in cross- 
sectional area of the specimens was found to be small during these 
tests, stress calculations were based upon the initial area. Both 
longitudinal and transverse strains were recorded using the electrical 
strain gauges and instrumentation described previously. 

4.2.3 Relaxation Tests 

Relaxation tests were conducted using a Twing-Albert test 
machine with initial strain applied at a head rate of 20 in. /min. 
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Figure 4. Creep Testing Apparatus 
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(50.8 cm./min.). This apparatus is shown in Fig. 5. The strain level 
was found to remain constant for the complete duration of each test. 

The same load cell and strain gauge instrumentation as described for 
the creep tests was used. However, the signal from the load cell was 
amplified and recorded continuously on a Hewlett Packard (7100 B) dual 
channel strip chart. 

4.2.4 Loading-Unloading Tests 

The loading-unloading tests employed essentially the same equip- 
ment as the constant strain rate tests. However, a considerable amount 
of this data was recorded on cassette tape and transmitted to the 
IBM/370 facility at VPI&SU for processing using the CB^ data acquisi- 
tion system [26] which was under development at the time. The system 
utilized the signals from the Instron load cell and Vishay amplifiers. 
Thus the essential components of the experimental apparatus were 
identical to those used in the constant strain rate tests. 

4.2.5 Elevated Temperature Tests 

Relaxation tests at various temperatures above and below the 
glass transition temperature of the adhesives were conducted on the 
Instron Machine in combination with an Instron environmental chamber. 
The apparatus used is shown in Fig. 6. Special test fixtures were 
fabricated to position the specimen inside the chamber. A copper 
cooling coil was affixed to the top fixture to dissipate heat away 
from the load cell. The axial strain was recorded using an Instron 
high temperature extensometer (Model G51-14A) which was capable of 
50% (0.5 m/m) strain. It is shown mounted on a specimen, in the 
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Figure 5. Relaxation Testing Apparatus 
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environmental chamber, in Fiq. 7. The signals from the extensometer 
and the load cell were amplified and recorded as in the relaxation 


tests described above. 



Figure 7. Specimen Mounted in Environmental Chamber 
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Chapter 5 

^■CONSTANT STRAIN RATE BEHAVIOR 

This chapter is concerned with the experimental results from the 
constant strain rate tests and the ability of the constitutive equa- 
tions. reviewed in Chapters 2 and 3 to predict this behavior. The 
first two equations considered are the Ramberg-Osgood models, Eqs. 
(2.2.7) and (3.3.1), and their respective rate forms, Eqs. (2.2.10) 
and (3.1.2). The effect of the stress-whitening; phenomenon on these 
equations is also presented. The modified Bingham model, Eqs. (2.4.3), 
is considered next; and the effect of the stress-whitening phenomenon 
on this model is observed through the use-of the nonlinear model, 

Eqs. (3.4.6), 

5.1 Constant Strain Rate Results 

The stress-strain behavior of Metlbdnd 1113 is shown in Fig. 8 
for four different strain (Head) rates. Linear elastic behavior with 
little rate effect was found for low stress levels. However, the 
elastic limit was found to be rate dependent. Above the elastic limit, 
significant rate effects were observed which can be attributed to 
viscoplastic behavior. The maximum stress was also found to be sig- 
nificantly rate dependent. The failure strains, while rate dependent, 
did not follow a consistent pattern. Similar results were observed 
for Metlbond 1113-2 (Fig-. 9). It can be observed, by comparison, that 
the modulus, elastic limit stress, and maximum stress are lower for 
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= 0.366* 
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* AVERAGE VALUES 
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Figure 8 . Stress-Strain-Strain Rate Response of Metlbond 1113 
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Figure 9. Stress-Strain-Strain Rate Response of Metibond 1113-2 
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Met! bond 1113-2. This is attributed to the reinforcing effect of the 
synthetic carrier cloth present in Metlbond 1113. The strain rates 
recorded in Figs. 8 and 9 were calculated from the slope of the 
strain-time data and were found to be constant up to local yield 
{i.e., the initiation of crazing or microcracks), after which the 
strain rate decreased slightly due to a. small reinforcing effect of 
the strain gauge. 

Figs. 10(a) and (b) are photographs of specimen surfaces before • 
and after testing, respectively. The failed specimen is typical of 
all failures observed during this investigation. The failure surface 
is perpendicular to the- loading axis and has the general appearance of 
a brittle tensile failure. Essentially no permanent necking was 
observed. 

5.1.1 Stress -Whitening Phenomenon 

In the rate tests, a stress-whitening phenomenon was observed to 
occur at stresses and strains well below the maximum stress levels. 

This phenomenon- was similar to that observed by Shouldberg and Lang 

[30] during their investigation of polymeric materials. It is believed 
that the stress-whitening observed was a crazing phenomenon which is 
frequently observed in polymers [31]. Crazing in a material may be 
defined as the occurance of localized highly elongated regions (crazes) 
whose density may vary from zero (true crack) to that of the virgin 
material depending upon conditions during their formation [32]. Hull 

[31] states that crazes formed in a uniaxial tensile stress field have 
.a shape similar to a crack, and the plane of the craze is at right 
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a. Untested Specimen 


b. Failed Specimen 


Figure 10. Specimen Surfaces of Me tl bond 1113-2 
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angles to the stress axis. Hull also states that in a complex stress 
field, the plane of the craze is normal to the maximum principal 
tensile stress. Since complex stress fields often exist in areas of 
local yield, it would then be possible for a specimen in which areas 
of localized yield occur during a uniaxial tension test to have crazes 
both normal and at various angles to the loading direction. 

The stress-whitening may be seen in the photograph shown in 
Fig. lOCb) for a specimen of Metlbond 1113-2. An untested specimen 
was shown in Fig. 10(a) for comparison. The stress-whitening pattern 
observed in Fig. 10(b) is very similar to a photograph presented by 
Hull [31] showing the crazes which formed prior to the maximum stress 
during a tensile test on polystyrene. Fig. 11 is a photograph of a 
failed specimen surface at approximately 30X. it can be observed that 
the crazes appear both normal and at various angles (vertical markings 
are from sanding) to the loading direction. These angles range from 
approximately ± 45" to ± 90° with respect to the loading axis. It is 
interesting to note that a considerable amount of these angles are 
similar to the angle (+ 54.1° to the- loading axis) at which slip 
bands occur in many metals in uniaxial tension [8]. ■ 

The time, stress and strain levels where the whitening occured • 
in the rate tests can be found in Table 1. The stress level at the 
occurence of stress-whitening was rate dependent, but the strain level 
was nearly constant. It is believed that the observed stress- 
whitening effects can be attributed to crazing and that this process 
is evidence that a local damage or failure mechanism occurs well in 
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Figure 11. Failed Specimen Surface at 30X 
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Table 1. CONSTANT STRAIN RATE PROPERTIES. 


Specimen 

Symbol s 
for 
Figure 
13 

Strain 

Rate • 

• 

e 

(m/m/sec) 

Elastic 

Modulus 

E 

(MPa) 

Elastic 

Limit 

Stress 

e 

(MPa) 

Elastic 

Limit 

Strain 

^0 

(m/m) 

B-4-3 

13(a), O 

6.81x10"® 

‘2161 

22.48 

0.0104 
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advance of gross fracture. As such, it is an important characteristic 
of this material. 

For a material in which the density is lower in the crazed area, 
the material in this location has undergone a type of phase change.' 

To' determine whether the crazed areas observed in the adhesives were 
true cracks (zero density) or material of less density, the microscopic 
picture shown in Fig. 11 was taken with the light source at'anangle 
to the specimen surface. The photograph indicates that no shaded 
valleys (true cracks) are visible. This would indicate' that the 
crazed areas of the adhesives are not true cracks, but areas of lower 
density.- Although subsurface studies were not conducted, it is be- 
lieved that this condition exists throughout the specimen since the 
primary stress field (uniaxial stress) responsible for this phenomenon 
is constant across the thickness. This change of density (or phase) 
is shown to influence the adhesives ' constitutive equations. 

Table 1 also gives further information regarding the properties 
found for the adhesives in the rate tests. The elastic moduli, 
Poisson's ratios for the initial elastic region, elastic limit stresses 
and strains, and maximum stresses may be found therein. 

Variations in Poisson's ratio with strain rate and with time are 
shown in Fig., 12. The results are shown in terms of a nondimensional 
time which is. based upon, the time required for failure in each test. 

The values for the failure times are given in Table 1. The region at 
which stress-whitening occured during the rate tests is also shown. 

It is believed that the occurence of stress-whitening is the reason 
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a. METLBOND 1113 


(SEE TABLE 1 FOR SYMBOL IDENTIFICATION) 



Figure 12. Poisson's Ratio During Constant Strain Rate Tests 
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for the general decrease in Poisson's ratios' with time. This decrease 
is Poisson's ratio is much. less pronounced in Metlbond 1113 as 
compared to Metlbond 11 13-2. This is attributed to the stabilizing 
effect of the carrier cloth present in the former adhesive. 

5.2 Ludwik's Equation 

Before investigating the ability of the constitutive equations 
to model .the constant strain rate behavior, certain essential 
empirical relationships describing the rate behavior of the elastic 
limit stress and strain and maximum stress 'observed in the adhesives 
need to be discussed. 

The elastic limit stress, 0, and elastic limit strain, e^, have 
been defined as the value of stress and- strain separating the elastic 
region from the viscoplastic region. These values were taken from 
data such as that shown in Figs. 8 and 9 to be the upper limit of. 
linearity for a given rate. Values for 0 and eq are recorded in 
Table 1. The maximum stress, Y, has been defined to be value of stress 
at the initiation of perfectly plastic flow. While these definitions 
are to some extent arbitrary, they do provide for consistent compari- 
son of experimental results with an empirical relationship proposed 
strictly for yield behavior by Ludwik [27]. 

As proposed by Ludwik, the variation in yield stress, Oy, with 
strain rate may be written as, 

(Jy = Oy' + CTy"log f, (5.2.1) 

where oy', Oy", and e' are material constants. 
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The elastic limit stress and yield-stress are identical. It is 
proposed,- therefore, that the variation of elastic limit stress and 
strain, and maximum stress with strain rate follow similar expressions. 
These can be written as. 


0 = 0'+ e"log l-r 

(5.2.2) 

£ 


0 = Co' + ^o"''og |T- 

(5.2.3) 

Y = r + Y"log fr 

(5.2.4) 


e 


where additional material constants are defined accordingly. 

Fig. 13 shows a comparison of experimental results for e and eq 
as a function of e compared with Eqs. (5.2.2) and’ (5.2.3), respective- 
ly. It may be observed that the equations fit the experimental data 
quite well and could be used to interpolate the data to other strain 
rates within the range of data. The variation of the maximum stress 
with strain rate is shown in Fig. 14. These results compare very well 
with Eq. (5.2.4). 

5.3 Modified Ramberg-Osgood Models 
5.3.1 RAMOD-1 

The first Ramberg-Osgood model studied is Eq. (2.2.7) which is 
repeated below as, 

e = + Ep = ^ + (5.3.1) 

and is called RAMOD-1 (Ramberg-Osgood. Model - 1). The modulus, E, is 
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b METLBOND 1113-2 


Figure 13. Variation of Elastic Limit Stress and Strain and 
Comparison with Ludwik's Equation 
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a.METL80ND 1113 
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b METLBONO 1113-2 


Figure 14. Variation of Maximum Stress and Comparison with 
Ludwik's Equation 
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taken directly from the stress-strain curve. Values. for K and n are 
determined from a log-log plot of a versus ep for each strain rate. 

Using values of E, K, and n from each rate test, RAMOD-1 was 
used to describe the behavior of Metlbohd 1113 and 1113-2. This 
model is called RAMOD-T-A. From the data presented in Figs'. 8 and 9 
log-log plots of a versus sp were determined and fitted with a least 
squares approximation. These curves are presented in Figs. 15. It 
can be noted that these curves deviate from the linear' approximation 
at higher values of stress and plastic strain. 

Values for E, K, and n for each strain rate were determined as 
described earlier and are recorded in.Table 2. Figs. 16 and 17 are 
comparisons of RAMOD-1 -A to experimental, data for Metlbond 1113 and 
1113-2, respectively. Ludwik's equation- (Eq. 5.2.4) was used to 
determine the maximum stress. It can be -observed that the data for 
both adhesives is approximated reasonably well up to the maximum 
stress with this model. 

In the rate form of RAMOD-1, Eq. (2.2.10), the coefficient n 
was considered to be invariant with strain rate. Using an average 
value of n, n, from the rate data, Eq. (2.1.14) can be written as, 

t 

E = ~T^ + as^a*^ (5.3.2) 

and is called RAMOD-l-B. The coefficients a and b were determined 

✓ 

from a log-log plot of K versus e; coefficients c and d were deter- 
mined from a log-log plot of E, versus e. Values for these coefficients 
and n are recorded in Table 2. 
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Table 2. COEFFICIENTS FOR RAMBER6-0SG00D MODEL - 1 (RAMOD-1) COEFFICIENT. 


Strain Rate 
(m/m/sec) 

K 

(MPa“") 

n 

Kl 

(MPa""'') ' 

"1 

K2 

(MPa""2) 

"2 

a(sec'^MPa"*^) 

3.767^10“^^ 

b 

-1.779 

c(sec*^MPa) 

2388.1 

6.81x10"^ 

5.505x10-43 

6.222 

3.048x10"^^ 

5.734 

7.197x10"®^ 

15.816 

•d 

0.0077 

6.83x10~5 

3.251x10-'* 5 

7.402 

1.336x10"'' 4 

7.011 

1.829x10-25 

13.506 

■ f 

11, .34 

,6.55x10-4 

1.787x10-'15 

7.388 

7.374x10-15 

7.003 

5.629x10"^ 

11.118 

g 

0.0496 

6.68x10“^ 

. - 1 

8.d63xl0-''9 

1 9.087- 

1.368x10"''® 

8.967 

il. 962x10-25 

13.261 

F 

7.525 
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K 

(MPa"*^) 

n 

Kl 

(MPa’"l) 


K2 

(MPa‘*^2) 

n2 

a(sec*^MPa-*^) 

5.257x10“!® 

Strain Rate 
(m/m/sec) 

"1 

b 

-1.606 




c(sec‘lMPa) 

2069.6 

7.00x10-5 

2.284x10- ‘ « 

5.574 

5.628x10-' • 

5.300 

1 .000x10-22 

12.697 

d 

0.004 

6.75x10"® 

2.188x10-11 

5.434 

1.080x10-1® 

4.962 

7.943x10-2® 

16.183 

'f 

9.91 

7.05x10-4 

1.856x10-1" 

6.537 

5.099x10-1® 

6.245 

?.264x10"®2 

18.000 

•• g 

0.0539 

7.70x10"® 

4.405x10-1^ 

7.959 

1.033x10-1® 

7.723 

5.1 96x1 0"^4 

12.686 

F 

6.376 


b. Metlbond 1113-2 
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Figure 16.. Stress-Strain-Strain Rate Behavior of Metlbond 1113 
and Comparison to RAMOD-V-A 
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Figure 17. Stress-Strain-Strain Rate Behavior of Metlbond 1113-2 
and C(OTparison to RAMOD-l-A 
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RAMOD-l-B did not model the rate' behavior of the adhesives very 
well and the theoretical curves are not shown for this reason. The 
main difficulty' is that rr is not invariant with strain rate. Thus the 
average value of n predicts very large plastic strains for the lowest 
strain rate and small plastic strains for the highest strain rate. 

This, of course, does not agree with the experimental data for the 
range of strain rates investigated. 

It is proposed to modify the form of RAMOD-l-B to allow for the 
rate dependency in n. The proposed form is, 

- ^ j. ,‘b fe^ /c o 

e - ^TT + ae a (5.3.3) 

ce^ 

and is called RAMOD-l-C. The coefficients a, b, c, and d are identical 
to the previous model; the coefficients f and g (also recorded in 
Table 2) were determined from a log-log plot of n versus * using a 
least squares approximation. Although this model was a definite im- 
provement over RAMOD-l-B, it did not adequately model the behavior of 
the adhesives for the complete range of strain rates. The main 
deterent here is that the proposed form for the rate dependency in n, 

n(e) = fe^ (5.3.4) 

does not adequately model the experimental data. 

5.3.2 Bilinear Form of RAMOD-1 

Re-examination of Figs. 15(a) and (b) indicates that the results 
presented in these figures may be better represented by two straight 
lines for each strain rate. Figs. 18(a) and (b) are the bilinear 



log CT {G’in MPa) 
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approximations. The stress value at which these Tines intersect has 
been defined as a*. Thus, RAMOD-1 has the form, 

a 

e = ^+K-|a ' 0<a<a* 

(5.3.5) 

' CT "2. 

e = ^ + ^20 0* < a < Y 

and appropriately called RAMOD-l-BL. .Values for the bilinear co- . 
efficients for K and n (K-j, K2,- n^ and n 2 ) are recorded in Table 2. 
Results using RAMOD-l-BL are compared to the experimental data in 
Figs. 19 and 20 for Metlbond 1113 and 1113-2, respectively. This 
•mode.l is an improvement over RAMOD-1 -A in which values of K and n 
were used. These improvements are a better representation of the slope 
throughout and the prediction of larger strains corresponding to the 
values of stress prior to the maximum stress. 


5.3.3 RAMOD-2 

A Ramberg-Osgood type model was proposed, Eq. (3.1.1), in which 
the plastic strain, ep, is a function of the over-stress above the 
elastic limit stress. This equation has the form. 


CT 

^ ~ £ 0 < a < Q 

e .= ^ + K(a - g)*^ 6 < 0 < Y 


(5.3.6) 


and is called RAMOD-2. Values for K and n are determined from a log- 
log plot of (0 - e) versus ep for each strain rate. As in the previous 
model, variations in E with strain rate denote viscoelastic effects. 
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Figure 19. Stress-Strain-Strain Rate Behavior of Metlbond 1113 
and Comparison to RAMOD-l-BL 
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Figure 20. Stress-Strain-Strain Rate Behavior of Metlbond 1113-2 
and Comparison to RAMOD-l-BL 
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while the parameters K and n represent variations in work-hardening 
characteristics with strain rate. 

Proceeding in' the same manner as RAMOD-1 , the first model in- 
vestigated was RAM0D-2-A which is actually written in the form of 
Eq. (5.3.6) and employs the values of E, K, and n for a given strain 
rate. Log-log plots of (a - 0 ) versus ep are shown in Fig. 21. A 
least squares approximation for these results yielded the values for 
K and n for each strain rate whichare recorded in Table 3. It can' 
also be observed in Fig. 21 that these results may be represented by 
a bilinear approximation. 

Results using RAM0D-2-A are compared to the experimental data 
in Figs'. 22 and 23. Ludwik's equations, Eqs. (5.2.2) and (5.2.4), . 

were used to determine the theoretical values for the elastic limit 
and maximum stresses, respectively. The results presented in Figs. 

22 and 23 indicate that this model predicts the behavior of both ad- 
hesives reasonably well up to the maximum stress. However, it would 
be desirable to have a better representation of the slope throughout 
and the prediction' of more plastic flow at the higher stresses. 

In the derivation of the rate form of RAMOD-2,' Eqs. (3.1.2), 
the coefficient n was also considered to be invariant with strain 
rate. Using an average value of n, rT, Eqs. (3.1.2) can be written as, 

e = 0 < o <'e 

ce^’ 

(5.3.7) 

0 . . b / . 

-"TT + ae (ff - e) 0 < 0 < Y 
ce 
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Table 3 . COEFFICIENTS FOR RAMBERG-OSGOOD MODEL - 2 (RAMOD-2) COEFFICIENT.- 


Strain Rate 
(iti/m/sec) 

■ K 

(MPa"") 

n 

Kl 

(MPa“'^1) 

"1 

•^2 „ 
"■rio 
(MPa 2) 

"2 

a(sec®^MPa“") 

2.778x10"5 

b 

0.102 

c(sec'^MPa) 

2388.1 

6.81x10-6 

7.318x10-6 

2.316 

1.370x10^^ 

2.045 

1.715x10“''3 

7.958 

d, 

0.0077 

6.83x10-6 

9.428x10-6 

2.247 

1. 437x1 0"5^ 

2.053 

2.457x10"^^ 

7.210 

■ • f 

• 1:885 

6.55x10-4 

2.370x10"^ 

1.960 

3.818x10“^ 

1.733 

3.114x10"^ 

4.912 ■ 

g ■ 

-0.0162 

6.68x10-3 

1.175x10-5 

2.141 ■ 

1. 689x1 0"5 

■ 1.972 

8.647x10-9 

4.574 

rT * 

2.166 
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Strain Rate 
(m/m/sec) 

K 

(MPa-") 

' n 

Kl 

(MPa“"l) 

"1 

K2 . 

(MPa'"2) 

n2 

a(sec6MPa ") 

1^430x1 0-4 

b 

0.233 

c(sec*^MPa) 

2069.6 

7.00x10"^ 

4.961x10"^ 

2.568 

7.365x10-6 

2.400 

1.073x10"’'^ 

8.969 

d 

o;oo4 

6.75x10"^ 

2.968x10"^ 

1 .973 

4.544x10-5 

1.772 

1.505x10-^^ 

10.419 

f 

1.500 

7.05x10-4 

4.115x10-5 

,1.824 

5.260x10-5 

1.689 

4. 062x10"’^ 3 

12.478 

g 

-0.0379 

7.70x10-3 

2.779x10-5 

1.957 

3.397x10-5 

1.848 

2.971x10-''3 

8.227 

rT 

2.081 


b. Metlbond 1113-2 
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Figure 22. Stress -Strain-Strain Rate Behavior of Met! bond 1113 
and Comparison to RAM0D-2-A 
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Figure 23. Stress-StrainrStrajn Rate Behavior of Metlbond 1113-2 
and Comparison to RAM0D-2-A ^ 
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and termed RAM0D-2-B. The coefficients a and b were determined from 
a log-log plot of K versus £ as shown in Fig. 24. The coefficients c 
and d were determined from a log-log plot of E versus e (Fig. 24). 
Values for these coefficients and n are recorded. in Table 3. 

The theoretical results using RAMOD-l-B did not adequately 
represent the rate behavior for either adhesive. This is attributed 
to the fact that although n is nearly a constant, RAMOD-2 is. very 
sensitive to the small variations in n with strain rate. 

It is proposed to modify RAM0D-2-B to account for the rate 
dependency in n. The proposed form is, 

a ' - 

e = —TT 0 < a, < 0 • 

. ce . 

(5.3.8) 

e =“ + a£^(<7 - 0)^^ ' 9 '< a < Y’ 

Ce“ 

and titled RAM0D-2-C. The coefficients f and g were calculated from 
log-log plots of n versus e i.n'Fig. 25 and are recorded in Table 3. 
Figs. 26 and 27 are comparisons of RAM0D-2-C with the experimental 
data for Metlbond 1113 and 1113-2, respectively. This model was a 
definite improvement over the previous rate form, RAM0D-2-B, and 
actually predicts the rate behavior better over all than RAM0D-2-A. 

From a programming aspect RAM0D-2-C is ideal in that only one 
parameter, e, needs to vary in the input data in order to study rate 
effects . 

Models RAM0D-2-A and RAM0D-2-C are adequate in predicting the, 
strain rate response of the adhesives. However, there is still a need 
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Figure 24. Log K and Log E Versus Log e for RAM0D-2-B and 
RAM0D-2-C 
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Figure 26. Stress-Strain-Strain Rate Behavior of Met! bond 1113 
and Comparison to RAM0D-2-C 
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Figure 27. Stress-Strain-Strain Rate Behavior of Metlbond. 1113-2 
and Comparison to RAM0D-2-C 
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for an improvement' both the representation of slope and the prediction 
of larger plastic strains correspond to stresses near the maximum 
stress. This suggests -the bilinear approach used in RAMOD-1. 

5.3.4 Bilinear Form of RAMOD-2 

Fig. 28 is the bilinear approximations of log (a - e) versus 
log £p. The stress magnitude at the intersection of the two straight 
lines for each strain rate has again been defined as a*. The bilinear 
values of K and n generated for each strain rate are given in Table 3. 
Thus RAM0D-2-BL has the form, 

e = |- 0<a<6 

.£ = 1+ K-|{a - e)^ e < a. < a* (5.3.9) 

_ np 

£ = g-+K2{a-0) <t*<0<Y 

Results using RAM0D-2-BL are compared to the experimental data 
in Figs. 29 and 30 for Metlbond 1113 and 1113-2, respectively. Obser- 
vation of these results indicate that a better representation of slope 
throughout and the prediction of larger plastic strains corresponding 
to stresses just prior to the maximum stress are realized. 

5.3.5 Stress-Whitening-Stress and a* 

In the bilinear models of both RAMOD-1 and RAMOD-2, o* has been 
defined as the stress at the intersection of the bilinear curves. It 
is interesting to compare these values to each other, and also to the 
value of the stress -whitening stress observed experimentally. These 
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Figure 29, Stress-Strain-Strain Rate Behavior of Metibond 1113 
and Comparison to RAM0D-2-BL 
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STRAIN (m/m) 

Figure 30. Stress-Strain-Strain Rate Behavior of Metlbond 1113-2 
and Comparison to RAM0D-2-BL 
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results are presented 1n Table 4. It can be observed that the a* 
values for both stress and strain are very close for the two models. 
These values are also slightly lower than the stress-whitening values. 
It should be noted that the point at which stress-whitening occured 
during the rate tests was recorded by visual observation. Thus, the 
true point of initiation of the whitening may have occured earlier and 
was not readily visible. The strain values corresponding to a* were 
relatively constant for all rates. It is apparent that in the region 
of stress-whitening, a change of phase occurs both physically in the 
material and analytically in the bilinear models. 

5.4 Modified Bingham Model 

The constitutive equations for the modified Bingham model (Eqs. 
2.4.1), which is shown in Fig. 1(a), are 

• .0 

e = "^ cr < 0 

(5.4.1) 



A stress such that a > Y is not allowed. For a constant strain-rate 
test, the solution of Eqs, (5.4.1) is Eqs. (2.4.3). This solution can 
be written as, 

o(e) = Ee a ^ 0 

(5.4.2) ^ 

a(s) = 0 + ERt(1 - e 6 < a < Y ’ 

where t = y/E is the relaxation time. 



Table 4. o* AND STRESS-WHITENING STRESS. 


Strain Rate 

Stress Whitening 

RAMOD- 

1-BL 

RAM0D-, 

2-BL 

Stress 

Strain 

a* 

Strain 

0* 

Strain 

(m/m/sec) 

(MPa) 

(m/m) 

(MPa) 

(m/m) 

(MPa) 

(m/m) 

6.81 X 10-6 

45.57 

0.0330 

44.58 

0.0293 

44.11 ■ 

0.0278 


49.64 

0.0332 


0.0295 

47.83 

0.0284 


53.37 

0.0348 

51.84 

■MEiiH 

51.44 



1 

[ 

I 

1 

1 

— 

— 

55.44 

0.0295 


a. Met! bond 1113 


Strain Rate 
(m/m/sec) 

Stress Whitening 
Stress Strain 
(MPa) (m/m) 

RAMOD-l-BL 
a* Strain 

(MPa) (m/m) 

RAM0D-2-BL 
a* Strain 

(MPa) (m/m) 

7.00 X 10-6 

39.23 

0.0378' 

38.77 

0,0345 

38.03 

0.0318 

6.75 X 10"^ 

42.47 

0.0377 

41.30 

0.0318 

41.31 

0.0308 


45.37 

0.0368 

mmm 

0.0321 

44.68 

0.0310 

7.70 X 10-^ 

1 

I 

1 

1 

1 

46.97 

0.0314 

47.08 

0.0304 


b. Met! bond 1113-2 
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Eqs. (5.4.2) were fitted to the rate dependent stress-strain 
behavior presented in Figs. 8 and 9. These results are shown in Figs. 
31 and 32 for Metlbond 1113 and 1113-2, respectively. Ludwik's 
equations for the e and Eq (Eqs. (5.2.2).and (5.2.3)) and Y (Eq. 
(5.2.4)) were used to determine the theoretical values for the elastic 
limit and maximum stress, respectively. With this procedure, it was 
possible to accurately represent the rate dependent behavior with the 
modified Bingham model. An accurate representation of slope is 
observed throughout. This model also allows perfectly plastic behavior 
after the maximum stress is reached which is in agreement with the 
experimental results. 

5.4.1 Rate Dependent Viscosity Coefficient 

It was found that a single relaxation time for each adhesive 
was inadequate to fit full range of the strain rate results. There-’ 
fore, a relaxation time (or viscosity coefficient, since u = tE) was 
determined to give the best fit at a given rate; The relaxation times 
necessary to achieve the close approximations obtained in Figs. 31 and 
Figs. 32 are plotted in Fig. 33. It is interesting to note that a 
single linear curve can approximate the data for both adhesives when 
plotted on log-log scales. The values for the relaxation times as 
well as the viscosity coefficients for each strain rate are presented 
in Table 5. It can be observed that as the strain rate increases an 
order of magnitude the viscosity coefficient decreases by an order of 
magnitude. 
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Figures!. Stress-Strain-Strain Rate Behavior of Metlbond 1113 and 
Comparison to Modified Bingham Model 
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Figure 32. Stress-Strain-Strain Rate Behavior of Metlbond 1113-2 
and Comparison to Modified Bingham Model 
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Figure 33. Relaxation Time Versus Strain Rate forfModified 
Bingham Model 


Table 5. RATE DEPENDENT RELAXATION TIMES 
AND VISCOSITY COEFFICIENTS. 


, & 

T 

y 

(m/m/ sec) 

(sec) 

(MPa - sec)- 


1800.00 



, 173.85 


6.55 X 10-4 

18.15 

4.07 X 1q4 

6.68 X 10-3 

1.85 


a 

. Metlbond 1113 

• 

e 

T 


{ m/m/ sec) 

(sec) 


7.00 X 10-6^ 

1816.80 

■ 3.67 X 105 

6.75 X 10-5 

180.20 

3.64 X 1Q5 


16.49 


7.70- X 10-3 

1.45 

2.93 X 103 


b. Metlbond 1113-2 
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The constitutive form of the modified Bingham model, Eqs. ■ 
(2.4.3), does not specifically indicate, that v may be rate dependent. 
However, the analogous form of these equations, Eqs. ,(2.3.15), 
derived from incremental theory indicated this possibility. Other ■ 
investigators, such as Chase and Goldsmith [22], using mechanical 
models to predict the behavior of rate-sensitive plastic materials 
observed. this same phenomenon. The rate dependency in y is therefore 
considered to be an inherent property of rate-sensitive materials. 

5.4.2 Nonlinear Model 

In an attempt to mathematically describe the rate effect ob- 
served in the viscosity coefficient for the modified Bingham model, 
the nonlinear perturbation technique proposed by Davis [4] (Section 
3.2) was investigated. This technique did not show the viscosity co- 
efficient to be dependent on the total strain rate e, since the non- 
linear hypothesis, Eqs. (3.2.1), states that variations in y are 
dependent on the strain rate in the dashpot, ey. However, the non- 
linear model did yield some interesting information relative to the 
viscosity coefficient of the adhesives above and below the stress- 
whitening point. 

Postulating that the nonlinear effects were present in the 
dashpot only, the nonlinear constitutive equations for the modified 
Bingham model (Eqs. (3.4.6)) can be written as, 
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o(e) = Eq£ 


a < 0' 


0(e,Yp) = e + EoRt( 1 - + Yy(2EoRx(e - Eq) 


+ EoR^T^d 


-2(e ^ £o)/Rt 

e ■ ) 


e < a < Y . 


where t =• yo/^o is the linear relaxation time as before.' For 6 < ct < Y 
the second of Eqs. (5.4.3) can be written as ' 

<?{e>Yp) = cfQ YyS-| e < a < Y (5.4,4) 

where 

00 = 0 + EoRt( 1 - " ^o)/R^) (5 4^5) 

is the linear solution as before and 

. n o “2(e - £n)/RT 

S-| = 2EoRt(£ - £q) + EqR t (1 “ e ) (5.4.6) 

The term 0 (e,Yp) is considered to be the experimental data, and 
the coefficient Yy is allowed to vary to fit the experimental data 
exactly. Then Eq. (5.4.4) can be written as, 

<°exp -^°th6o> , (5 4 7 ) 

where c^thgo " ‘^o represents the theoretical curve from the linear 
solution. As an example of the implications of Eq.- (5,4.7), the 
parameters in this equation were calculated using the data for Metlbond 
1113-2 at a constant strain rate of e = 7.05 x 10“^ m/m/sec (Fig. 9). 
These results are presented in Fig. 34. 
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5.4.3 Significance of Stress-Whitening Stress 

It is interesting to note the point (in Fig. -33) at which the ‘ 
linear solution and the experimental results agree exactly (i.e.,- 
Yy =0). The stress and strain values, 43.0 MPa and 0.305 m/m, 
respectively, are slightly lower but very close to the stress- • 
whitening values for this rate (i.e., 45.37 MPa and 0.0368 m/m, 
respectively). Similar -results can be obtained for the other strain 
rates. It' can also be observed that the value of Yy changes sign in 
this transition (stress-whitening) region. Therefore, the perturba- - 
tion analysis shows that. the value of the viscosity coefficient prior 
to the stress-whitening region is lower than that predicted' by the 
linear solution and higher above this region. This indicates that co- 
efficients present in the modifie.d Bingham could vary above and below 
the stress -whitening region in order to obtain a more accurate • 
representation of the experimental results. This same phenomenon was ' 
observed in the modified Ramberg-Osgood models. 

Although this nonlinear analysis yielded interesting information 
relative to the linear solution above and below the stress -whitening 
region, it has not been considered further in this study since the 
nonlinear effects present are very small. That is, the linear solution 
of the modified Bingham model adequately predicts the rate behavior. 

5.5 Models in Advanced laminate Analysis 

A goal of this study has been to develop analytical models 
capable of predicting the strain rate behavior of Metlbond 1113 and 
1113-2. These models should also be easily adaptable for use in 
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advanced analysis of material systems. Recent investigators [28, 

29, 30] have used various modified Ramberg-Osgood equations in. 
advanced analysis of composite laminates. Hashin, et al. [28] use a 
modified form of RAMOD-1 , written as,: 



(5.5.1) 


where ay is a curve fitting parameter, to describe the behavior of the 
matrix material in boron/epoxy laminates. Kiblef [29] also used the 
form of Eq. (5.5.1) to model the nonlinear behavior of the transverse 
extensional and inplane shear stress^strain relationships for uni- 
directional, fiber reinforced laminae. In an effort to describe these 
same latter two relationships, G. Renieri [30] is currently using a 
model in the form of RAMOD-1. The model is being used in an incre- 
mental loading scheme in conjunction with the finite element method 
in order to predict the stress distributions in composite laminates. 
RAMOD-1 was considered since it gave an accurate representation of the 
slope which is essential in the incremental analysis. For example, 
using RAMOD-1, the slope would be. 


dg _ 
de 


E ■ 

1 + EnKo""^ 


(5.5.2)- 


and similar equations would prevail for the various forms' of RAMOD-1 
and RAMOb-2. For stresses below the elastic limit, the slope, da/de, 
given by Eq. (5.5.2) is essentially the .modulus E. Above the elastic 
limit the slope begins to decay exponentially from a value of E. 
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Although the literature has not produced such evidence, the 
modified Bingham model could also be used in advanced laminate analysis. 
This model would be ideally suited for the analysis of material systems 
using Metlbond 1113 or 1113-2. For incremental formulations the 
slope using the modified Bingham model would be. 


d£ ■ 


(5.5.3) 


which shows that the slope is E at the elastic limit (e = Sq) and 
the n decays exponentially, as expected, for higher strains. 

The modified Bingham model could also be used to model a single 
nonlinear stress-strain curve (i.e., rate data not available). For 
example, Eqs. (5.4.2) can be written as. 


CT = Ee 


a < 8 


a = 6 + EKod 





e < a < Y' 


(5.5.4) 


where Kq is a maternal constant determined empirically to give the 
‘best fit to the experimental data. Equations of this form may be 
useful in modeling the nonlinear behavior of composite laminate in 
future -investigations. 



Chapter 6 

CREEP, RELAXATION, AND UNLOADING BEHAVIOR 

This chapter is concerned v/ith the ability of the modified 
Bingham model to predict the behavior of Met] bond 1113 and 1113-2 
for loading conditions other than constant strain rate. Creep results- 
are presented initially. The relaxation results are then reviewed. 

The final section of this chapter presents the experimental and 
theoretical results for two, more complex strain histories. 

6.1 Creep Results 

Creep results are shown in Figs. 35 and 36 for Metlbond 1113 and 
1113-2, respectively. For these tests the levels of stress were 
purposely larger than the elastic limit values found in the rate tests 
as little time effects were expected below these levels. 

It may be seen from Figs. 35 and 36 that the adhesives exhibit 
a delayed failure phenomenon. That is, in the creep tests, failure was 
observed to occur after a period of time at a particular stress level. 
The time to failure decreased with increasing stress level; these times 
are recorded in Table 6. These results are similar to the delayed 
failure phenomenon found for polycarbonate by Brinson [23]. Delayed 
stress-whitening was also found to occur and is evidence that both 
local and gross failure are time dependent. 

Variations in Poisson's ratio with stress level and time during 
the tests for both materials are shown in Fig. 37. These results are 
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Figure 36. Creep Behavior of Met! bond 1113-2 
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Table 6. CREEP PROPERTIES. 


Adhesive 


Specimen 


Symbol s 
for 

Figure 37 


Stress 

(MPa) 


Creep to 
Failure 
Time, tp 
(min) 


Metlbond 

1113 


B-6-2 

B-6-5 

B-5-3 

B-5-4 

B-5-5 


37(a), 

37(a), 

37(a), 

37(a), 


a 

□ 


o 


40.64. 
42.79 
47.22 
51.02 
56. 4¥ 


410 

260 

26.75 

3.35 

0.046 


Metl bond 
1113-2 


A-6-3 

A-6-1 

A-6-2 


37(b), □ 
37(b), 9 
37(b), O 


39.09 

41.02 

42.54 


54 

14 

5.9 

































0 0.2 OA 0.6 0.8 l.O 


t/tp 

a. METLBOND ill3 


(SEE TABLE 6 FOR SYMBOL IDENTIFICATION) . 



0 0.2 0.4 0.6 .0.8 1.0 


t/tp 

b. METLBOND 1113-2 

Figures?. Poisson's Ratio During Creep Tests 
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shown in terms of nondimensional time v/hich is based upon the time 
required for failure in each test. As in the rate tests, Poisson's 
rattos generally decreased with time with Metlbond 1113 being more 
stable due to the presence of the carrier cloth. After stress- 
whitening occurred in the rate tests, Poisson's ratio v/as observed to 
decrease more rapidly. In the creep tests stress-whitening occurred 
either during loading (that is, nearly instantaneously) or very early 
in time with respect to the failure time. Therefore, it is felt that 
this phenomenon is responsible for the lower initial values of Poisson's 
ratio in the creep tests as compared to the rate tests. 

The creep response of the modified Bingham model is given by 
Eq. (2.4.7), 

— 0 

e(t) = — t + ^ a>9 (6.1.1) 

Equation (6.1.1) did not adequately model- the creep behavior of the 
adhesives. The discrepancy here is that the modified Bingham model 
predicts a linear response for creep which obviously does not correspond 
to the experimental data. However, the modified Bingham model was 
found to accurately predict the delayed failure phenomenon when used 
in conjunction with an equation for maximum stress proposed by 
Crochet [31]. 

6.1.1 Creep to Failure Behavior 

A possible rational mathematical characterization of delayed 
failure phenomenon has been proposed by Crochet. He assumed that the 
maximum stress would increase for increasing strain rate and proposed 
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the equation, 


Y(t) = A + B exp(-Cx) (6.1.2) 

where A, B, and C are material constants and x is a time dependent 
material property given by, - 

X = (6.1.3) 

In Eq. (6.1.3) and refer to viscoelastic and elastic strains, 
respectively. 

In the work of Crochet, the material behavior was considered as 
viscoelastic/plastic. That is, below the maximum stress, material 
behavior was considered to be viscoelastic. Stresses above the maxi- 
mum stress were not allowed and perfectly plastic flow was postulated 
once the maximum stress was reached. 

For creep in uniaxial tension, Eq. (6.1.3) can be written as, 

X = C(en'^ - en^) + (£22^^ “ £22^)^ 

^^ 33 '^ “ (6.1.4) 

The second term in Eq. (6.1.1) represents the elastic behavior. Thus 
the term (ei/. - £-, 1 ^)^ in Eq. (6.1.4) becomes, 

( 6 . 1 . 5 ) 

If a constant Poisson's ratio is assumed,- the second and third terras 


of Eq. (6.1.4) can be written as. 
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/ V E»2 , V E\2 2 

^^22 “ ^22 ) “ (^33 ” ^33 ) ^ 


Oq - 0 


( 6 . 1 . 6 .) 


Thus Eq. (6.1.4) becomes, 


X = 




(6.1.7) 


and by using this result, Eq. (6.1.2) becomes, 

Y(t) = A + B exp[-C'(go -.e)t]. (6.1.8) 

where C = — ^ is a constant containing various material 
parameters. In a- creep test with a > 0, Y(t) and are identical. 
Hence ,' 


Y(t) = A + B exp{C'[Y(t) - e]t} 


(6.1.9) 


and.t represents the time for creep to failure to occur under constant 
stress. 

The material constants A, B, and C were determined for Hetlbond 
1113 and creep to failure times for different stress levels were calcu- 
lated according to Eq. (6.1.9). These results are shown in Fig.. 38, 
together with the experimental data. As may be observed, close 
correlation between measured and predicted values was found. 

6.2 Relaxation Results 

Relaxation results are shown in Figs. 39 and 40 for Met! bond 
1113 and 1113-2, respectively. For these tests, the levels of strain 
were purposely larger than the elastic limit values found in the rate 
tests since little time effects were expected below these levels. 
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Figure 38. Comparison Between Creep to Failure Data and Theory 
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Figure 40. Relaxation Behavior of Met! bond 1113-2 
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The relaxation response of the modified Bingham model is given 
by Eq. (2.4.11), 

. a(t) = 0'+ (sq E a>0 . (6.2.1) 

Eq. (6.2.1) states that, regardless of the strain level Fq, as t 
approaches infinity, o(t) asymtotically approaches e. However, it can 
be observed that the experimental data for each strain level does not 
approach the same value of e for large time. It can also be observed 
that the initial condition, 

a(t = 0) = Fq E (6.2.2) 

does not adequately represent the data. • 

In an attempt to use Eq. (6.2.1) to model the relaxation behavior, 
the initial condition was modified as, 

' a(t.= 0+)”=“dQ’ (6.2.3) 

where Cq' is the initial- experimental value of stress for a given 
strain level. Various schemes were then tried in order to fit the 
data. The modified Bingham model was unable to predict the relaxation 
behavior for the entire time span present in the experimental data. 
However, reasonable results were realized for short time (t < 2 min.) 
by using the value of the viscosity coefficient during loading^ and 

^The loading rate was 20 in/min (50.8 cm/min) which is an order 
of magnitude, above the highest head rate used' in the rate tests. Thus, 
the data in Fig. 33 was extrapolated to the a strain rate of e 
7.00x10"^ m/m/sec and a value for v determined. 
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adjusting the asymtotic value of 9. 

6.3 Loading-Unloading Results 

The loading-unloading tests were conducted on the most recent 
set of adhesive panels, to compare the properties- of these panels to 
the previous panels, constant strain” rate tests were conducted at 
0.2 in/min (0.508 cm/min). These results are shown iii Figs. 41 and 
42 for Met! bond 1113A and 1113-2A where the (A) represents the; second 
set of panels.. 

The properties of Metlbond 1113-2A are identical to those for 
Metlbond 1113-2 when tested at the same rate. Although' higher values 
in modulus, elastic limit and maximum stresses are observed in Metlbond 
1113A, these values are not as high as those for Metlbond 1113.. This 
has been attributed to the possibility that the properties of the 
carrier cloth varied between these two sets of panels. 

Constant strain rate, properties of Metlbond 11 13A and 1113-2A 
tested at a head rate of 0.2 in/min (0.508 cm/sec) are presented in 
Table 7. The values of the viscosity coefficient used to approximate 
the experimental results in Figs. 41 and 42 can also be found therein. 
The stress-whitening phenomenon was observed to occur in the same 
manner as in the previous panels-. 

6.3.1 Strain History- 1 

The first strain history considered is constant strain rate 
loading and unloading cycles followed by constant strain rate to 
failure. Fig. 43 is a schematic representation of one cycle, the 
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Figure 41. Stress-Strain Behavior of Metlbond 1TI13A and 
Comparison to Modified Bingham Model 
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Figure 42. Stress-Strain Behavior of Metlbond ni3-2A and 
Comparison to Modified Bingham Model 



1Q2 



TIME, t 

Figure 43. Strain History - 1 and Stress Response 
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assumed stress response is also shown. 

For loading-unloading considerations, the constitutive equations 
for the modified Bingham model, Eqs. (3.3;3)-, are, 

for a < e - 


e 



or a < 0- 


e 



for a > e 
and a > 0 


(6.3.1) 


That is, unloading occurs elastically and the elastic limit stress is 
non-work-hardenihg as shown in Fig. 42. The various times shown in 
Fig. 43 indicate the points in time where changes from elastic to 
viscoplastic behavior or from loading to unloading behavior occur. 

Eqs.- (6.3.1) are solved between these times for the stress .response. 
Values of stress and strain generated at the end of a given time seg- 
ment are used as initial conditions for the next segment. 

For the range 0-< t < tQ, the strain is given by, 

e(t) = Rt (6.3.2) 

and the stress response is, 

- . ' a(t) = ERt . (6.3.3)' ■ 

The time tg is determined from Eq. (6.3.3) as, 

. t„ = Ir . (6.3.4) 

. and Eq is determined from Eq. (6.3.2) as. 
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^0 ” ^^0 

In the interval tQ, < t < t].. 


(6.3.5) 


• e(t) =00+ R(t - to) 
and 

a(t) = 6 + yR(l - e'^^^ " 


(6.3.6) 

(6.3.7) 


The time t-| is known from the given load histor^^. The stress and 
strain af t-| are given by Eqs. (6,3.7) and (6.3.6) as, 


0 -| = e + nR(l - ■ ”^0^^^) 


(6,3.8). 


and 

.£•1 = 00 -*- R(t^ - tp) (6.3.9) 

respectively. 

For the range t-| < t < t 2 ,- the strain is, 


e(t) = - R(t - t^) 

and the stress response given by. 


(6.3.10) 


' cr(t) ='a-| - ER(t ~ t-, ) (6.3.11) 

The time t 2 is also known from the given load history. The stress and 
strain at t 2 are given by Eqs. (6.3.11) and (6.3.10) 


02 = 01 - ER(t2 - t-| ) (6.3.12) 

and 

02 = E] - R(t2 - t^) (6.3.13) 

respectively. 
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During the interval t 2 < t < t 2 ', the strain is given by, 

e(t) = £2 + - t2)‘ (6.3.14) 

and stress response is, • 

0(t) = 02 + ER(t - t2) (6.3.15) 

The time, t 2 'j at which the elastic limit stress is reached during the 
second loading can be determined from Eq. (6.3.15) as, 

d ^ Op ' 

^2* ~ ER (6.3.16) 

and the strain at t 2 ‘ is determined from Eq. (6.3.14) as,- 

£2' = £2 + " ^2^ (6.3.17) 

For the interval t 2 ' < t < ty, the strain can be written as, 

• £(t) = e^' + R(f^ t2') (6.3.18) 

and the stress response is, 

0(t) = 6 + yR(l - e'^^^ " (6.3.19) 

The strain at the maximum stress can be determined from Eq. (6,3.18) 
in the form, 

Ey = £2' + R(ty - t2') (6.3.20) 


where ty is given by Eq. (6.3.19) as. 



(6.3.21) 


lUb 




(e + yR) 
-yR 


.Above ty the material is assumed to flow perfectly plastically to 
failure. For a number of cycles, the analysis proceeds in a similar 
manner to that above. 

Metlbond 1113-2A and 1113A were tested for one and two cycles of 
constant strain rate loading and unloading, respectively, followed by - 
constant strain rate to failure. The material properties observed 
during these tests are presented in Table 7- The modulus during load- 
ing was found to be that of the initial elastic modulus. However, the 
modulus during unloading was found to be lower than the initial- 
elastic modulus. This phenomenon was observed for all unloading, tests. 
As such, it is an interesting property of these materials. The reason 
for this behavior is not clearly understood at this time. It should, 
however, be a topic, of future research. Values for the loading and 
unloading moduli are presented in Table 7. 

Using the material properties presented in Table 7 and equations 
in the form of Eqs. (6.3.2) to (6.3.21), a value for the viscosity 
coefficient (given in Table 7) for each adhesive was determined which 
gave the best fit to the experimental results. The experimental and 
theoretical results are compared in Figs. 44 through 49. It can be 
observed that an accurate representation of the experimental results 
was achieved. It can also be observed that during subsequent loadings 
a finite value of strain’ is predicted theoretically for the zero state 
of stress. This is in agreement with the experimental results, and 



Table 7. CONSTANT STRAIN RATE LOADING-UNLOADING PROPERTIES. 


Adhesive 

Speqimen 

. Strain 
Rate . 

ie| 

' 'avg 
(m/m/sec) 

Viscosity 

Co- 

efficient 

y 

(MPa-sec) 

Elastic 

•Limit 

Stress 

0 

‘(MPa) 

Elastic 
Limit 
■ Strain 
eo 

(m/m) ■ 

Elastic 

Modulus 

During 

Loading 

E 

(MPa)- 

Elastic 
Modulus 
During 
1st Un- 
loading 
El 

(MPa) 

Elastic 
Modul us 
During 
2nd Un- 
loading 

^2 

(MPa) 

Maximum 

Stress 

Y 

(MPa) 

Metlbond 

1113A 

B-7-2 


3.22x10^ 

27,5 

0.0131 




47.0 

B-7-5 


3.50x10^ 

27.5 

0.0131 

• 2100 


11991 

47.5 

Metlbond 

1113-2A 

A-8-2 

6.95xl0-'^ 


24.3 

0:0122 

bhi 

«... 


45.3 

A-8-4 



24.3 

0.0123 

1975 

wm 


45.5 
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Figure 44. Strain History - 1 for Metlbond 1113-2A and 
Comparison to Theoretical Approximation . 
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Figure 45. Theoretical and Experimental Stress-Time Response of 
Metlbond ni3-2A for Strain History 1 
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Figure 46. Theoretical and Experimental Stress-Strain Response 
of Metlbond 1113-2A for Strain History - 1 
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Figure 47. Strain History - } for Metlbond mSA and Comparison 
to Theoretical Approximation 
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, Figure 48. Theoretical and Experimental Stress-Time Response of 
Metlbond 1113A for Strain History - 1 
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Figure 49. Theoretical and Experimental Stress-Strain Response 
of Metlbond 1113A for Strain History - 1 
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physically this can be attributed to the residual plastic strain pro- 
duced during the previous loading cycle.. 

6.3.2 Strain History - 2 

The second strain history considered is constant strain rate - 
constant strain cycles followed by constant strain rate to failure. 

Fig. 50 is a schematic representation of one cycle. The assumed stress 
response is also shown-. The various times shown in Fig. 50 indicate • 
the points in time where changes from elastic tO' viscoplastic behavior 
or from constant strain rate loading to constant strain loading occur. 
During constant strain rate segements, Eqs,. (5.4.1) are solved accord- 
ing to the initial conditions generated at the end of the previous seg- 
ment. During relaxation (constant strain) segements, Eq. (2.4.5) is 
solved according to the appropriate initial condition'. 

For the range 0 < t < t-| , Eqs. (6.3,2) to (6.3.9) apply. In the 
interval t-j < t < t 2 » the strain is given by 

e(t) = e-j (6.3.22) 

and the stress response as, 

c(t) = e + (o, - (6.3.23) 

The time t 2 is known from the loading history. .The' values of- stress 
and strain at t 2 are given by Eqs. ■ (6.3.23) and (6.3,.22) as, 

. (72 = 6 + W - “ ■tl)/^' (6.3.24) 


and 



STRESS, cr 
STRAIN, € 


115 



TIME, t 


Figure 50, Strain History - 2 and Stress Response 
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£2 ■ 

respectively. 

During the interval t 2 < t < t^, the strain is. 


(6.3.25) 


e(t) = g 2 + R(t - t 2 ) (6.3.26) 

and the stress response is given, by, 

a(t) = 6 + yR + (c 2 - (8 + v>R))e ^ ^ (6.3.27) 

The strain, Sy, corresponding to the maximum stress is given by 
Eq. (6.3.26) as, 

£y = Gg + R(ty ~ t2) • (6.3.28) 

where ty is determined from-Eq. (6,3.27) as. 


y - (9 t TiR) 


ty “ *2 - E|_ct 2 - -(0 + vR). 


(6.3.29) 


Above ty the material is assumed to flow perfectly plasticly to failure 
For a number of cycles, the analysis proceeds in a similar manner to 
that above. 

Hetlbond 1113-2A and 1113A were tested ..for one cycle and two 
cycles of constant strain rate - constant strain, respectively., 
followed by constant strain rate to failure. The material properties 
observed during these tests are presented in Table 8. 

Using the material properties presented in Table 8 and equations 
in the form of Eqs. (6.3.2) to (6.3.9) and (6.3.22) to (6.3.29), a 
value for the viscosity coefficient (given in Table 8) for each 



Table 8. PROPERTIES FOR STRAIN HISTORY-2 


Adhesive 

Specimen 

Strain 

Rate- 

^avg 

(m/m/sec) 

Viscosity 

Co- 

efficient! 

y 

(MPa-sec) 

Elastic 

Limit 

Stress 

0 

(MPa) 

Elastic 
Limit 
• Strain 
^0 

(m/m) . 

Elastic 
Modulus 
■ E 
. (MPa) 

Maximum 

Stress 

Y 

(MPa) 

0 

During 
1st Re- 
laxation 
(MPa) 

0 

During 
2nd Re- 
laxation 
. (MPa) 

Metlbond 

.1113A 

B-7t6. 

6.80x10"^ 

.3.70x10^ 

27.5 

0.0126 

2180 

47.0 
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Metlbond 
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6.85x10'^ 


24.3 

0.0122 


46.3 

34.0 
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adhesive was determined which gave the best fit to the experimental 
results. It was found that the value of y determined during the 
constant strain rate loading can be used in the governing equations 
for the relaxation phases. The values of e during the relaxation- 
phase are then determined in order to give the best fit. As in the 
relaxation results, Section 6.2, -this procedure yields the best ap- 
proximation to the experimental results for short durations of time. 

The' experimental and theoretical results are compared in Figs. 

51 through 56. The values of e used in the governing equations during 
the relaxation phases are determined from the experimental results 
presented in Figs. 52 and 55. These values are taken to be the 
asymtotic value of stress approached if the relaxation curves are 
extrapolated beyond their time segment. The asymtotic values of 6 
are also recorded in Table 8. Since the relaxation times are short 
(t < 2 min) an accurate representation of the relaxation phases are 


represented using the above' procedure. It can also be observed that 
overall agreement between theory and experiment is achieved. 
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Figure 51. Strain History - 2 for Metlbond ni3“2A and 
Comparison to Theoretical Approximation . 
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Figure 52. Theoretical and Experimental Stress-Time Response of 
Met! bond 1113-2A for Strain History - 2 . 
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Figure 53. Theoretical and Experimental Stress-Strain Response 
of Metlbond ni3-2A for Strain History - 2 
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Figure 54. Strain History - 2 for Metlbond ni3A and Comparfson 
to Theoretical Approximation 
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Figure 55. Theoretical and Experimental Stress-Time Response 
of Metlbond lllSA for Strain History - 2 
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Figure 56. . Theoretical and Experimental Stress-Strain Response 
of Metlbond 1113A for Strain History - 2 




Chapter 7 

ELEVATED TEMPERATURE BEHAVIOR 

The mechanical properties of polymeric materials have a; very 
strong dependence upon temperature. At least three distinct regions 
of behavior can be observed for a cross-linked polymer when the 
temperature is varied [31]. These are the glassy region, the transi- 
tion region, and the rubbery plateau. The region which is of most 
interest in viscoelasticity is the transition zone, as time effects 
are especially pronounced in this region. Time effects are only ob- 
served over a long period of time (many decades of time) in the glassy 
range, whereas events occur so rapidly in the rubbery range that they 
are difficult to observe in an actual test. 

The regions above can be established by conducting relaxation 
tests at various temperatures and then plotting the relaxation modulus 
versus temperature for a given time. For such tests, high-modulus 
glassy- type behavior will be characteristic of the glassy region, while 
low-modulus rubbery-type behavior will be characteristic of the rubbery 
region. Results for Metlbond 1113A and 1113-2A are presented in Figs. 
57 and 58, respectively. It was found that both adhesives exhibited 
definite glassy and rubbery regions. It can be observed that these 
regions encompass the same temperature ranges for both adhesives. 

This indicates that the carrier cloth had little effect as to the 
temperature dependence of these regions. 
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TEMPERATURE i®C) 

Figures?. Variation in Relaxation Modulus with Temperatlire for 
Metlbond 1113A 
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The region where the modulus has the greatest dependence- upon 
temperature (maximum magnitude of slope), is the transition region. 
Within this region a corresponding temperature, or narrow temperature 
range, is designated as the glass transition temperature Tg. The 
glass transition temperature is located near the beginning of the 
transition zone [32].- The value of Tg is approximately 90® C for both 
adhesives. 

7 . 1 Time-Temperature Superposition 

A statement of the time-temperature superposition principle is 
given by the equation [31], 

E(t',To) = ^^E(t = a-j-t',T) (7.1.1) 

where p represents the density of the material at temperature T and aj 
is the "shift factor." This equation indicates that the relaxation ■ 
modulus at the reference temperature Tq and time t' can be obtained 
from the relaxation modulus measured at temperature T and time t. On 
a plot of modulus versus logarithmic time, this amounts to a vertical 
shift of the magnitude PqTo/pT and a horizontal shift of the magnitude 
log ay. Although time-temperature superposition is applicable to 
other viscoelastic response (creep, dynamic, etc.), it is usually 
illustrated with' stress relaxation, since most of the early development 
was done with stress relaxation [33]. 

The relaxation results were reduced by the factor Tq/T (the factor 
Pq/p was neglected as is frequently done in the literature [34]) and 
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are shown plotted in Figs. 59 and 60- for.Metlbond 1113A and lllS-ZA, 
respectively. The reduced curves were shifted laterally to obtain a 
master- curve for a reference temperature of Tq = Tg 90° C. It can 
be observed that this procedure yieids a smooth continuous curve 
which represents the stress relaxation at 90° C over thirteen decades 
of time. Multiplication by the appropriate value of aj then estab- 
lishes the master curve at any other temperature and can thus be used 
to predict response at that temperature over the entire time scale. 

7.1.1 WLF Equation 

For many polymeric materials it has' been found that if the 
polymer's glass transition temperature is chosen as the reference' 
temperature, the shift factors are given, to a good approximation, by 
the Williams-Landel -Ferry (WLF) equation [3.1],- 

-Ci(T-T) 

Log aj = C 2 + (T - Tg) (7.1.2) 

where C-] = 17.44 and C 2 = 51.6. Eq. (7.1.2) is only valid from Tg up 
to approximately Tg + 100° C since Eq. (7.1.1) cannot be proven to be 
valid for temperatures lower than Tg [31]. However, most researchers 
believe the time-temperature superposition principle is valid below 
the Tg— exactly how far below no one has ascertained [32]. For this 
reason, results for T < 70° C are not shifted in Figs. 59 and 60. 

The values of a-p used to superpose the various curves in Figs. 

59 and 60 are compared with the WLF equation in Figs. 61 and 62, 
respectively. Although the WLF equation does not adequately fit the 


y , -;'ni )l '( 


li ;; 





REDUCED RELAXATION MODULUS, E{t) (MPa) 


13Q 



log, ©TIME (min.) 


Figure 59. Master Curve of Relaxation Modulus for T© = 90® C 
. for Metlbond 1113A • 
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Figure 60. Master Curve of Relaxation Modulus for T© = 90® C 
for Met! bond m3-2A 
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Figures!. Comparison Between Experimental Data and WLF Equation 
for Met! bond !!13A 
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Figure 62. Comparison Between Experimental Data and WLF Equation 
for Met! bond m3-2A 
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experimental results, the shapes of these curves are similar, indicating 
that C-| in the WLF equation should not be taken as C-| = 17.44. If the 
constant C-| is adjusted for each adhesive (as given in Figs. 61 and 
62), the curve for the WLF equation is brought into reasonably close 
agreement to the experimental results. For t < 70° C and t > 120° C 
the experimental shift factors begin to diverge considerably from the 
WLF equation, indicating that the WLF equation is invalid outside of 
these temperatures. It should be pointed out, however, that it is 
possible for a shift factor to exist even though the WLF equation is 
not satisfied [32]. 



Chapter 8 
CONCLUSIONS 

The present investigation has been concerned with the stress- 
strain, strain rate, creep, relaxation, yield and failure properties 
of two structural adhesives used in composite applications, loading- 
unloading and elevated temperature behavior have also been studied. 

In summary of the achievements of this study, the following con- 
clusions can be made. • . ' 

(1) Bulk "form testing is a viable means of obtaining mechanical 
properties of- structural adhesives. 

(2) Me-tlbond 1113 and 1113-2 adhesives exhibit both rate and 
time dependent behavior. It was observed that the elastic 
limit stress and strain, and maximum stress are rate depen- 
dent; this behavior is accurately predicted using an 
equation proposed by Ludwik for yielding. The viscosity 
coefficient, y, is also rate dependent.- This phenomenon is 
indicated by the analogous constitutive form of the modified 
•Bingham model derived from incremental theory, Eqs. (2.3.5). 

(3) A stress'-whitening phenomenon which -is observed in the 
adhesives at high stress levels is attributed to crazing. 

The material in the crazed area is less dense than the 
virgin material. The properties of the adhesives are dif- 
ferent after stress-whitening with a resulting increase in 
plastic flow and decrease in Poisson's ratio. 
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(4) Higher; strength Is observed in Met! bond 1113 as compared to 
Met! bond 1113-2 due to the presence of the carrier cloth. 
This carrier cloth also tends to stabilize Metlbond 1113 

as failure strains during the rate tests are less random 
and Poisson's ratios above the stress-whitening level are 
relatively constant. 

(5) It is shown' that constitutive equations postulated by early 
investigators based on experimental results can be derived 
from existing theories and mechanical models. 

(6) The proposed bilinear forms of the Ramberg-Osgood model, 

Eqs. (5.3.5) and (5. 3'. 9) accurately predict the rate 
behavior of the adhesives. These models indicate that the 
material coefficients are different before and' after stress- 
whitening due to changes in the material properties'. The 
bilinear forms predict an accurate representation of the 
slope of the stress-strain curves which is essential when 
incremental formulations are used in advanced analyses. 

(7) The modified Bingham model, Eqs. (5.4.1), accurately pre- 
dicts the rate behavior of the adhesives. This model gives 
an accurate representation of the slope without having- to 
change material coefficients before and after the s tress - 
whitening point. .This model also has the advantage of pre- 
dicting perfectly plastic flow once the maximum stress is 
reached; such behavior is exhibited by the experimental 
results. It is felt, therefore, that of all models 
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investigated, the modified Bingham model is the preferred 
one for representing the rate behavior of the adhesives. 

(8) The nonlinear form of the modified Bingham model, Eqs. 
(5.4.3), indicates that the coefficients of the linear model 
may vary due to the change in material properties before and 
after the stress-whitening point. The nonlinearities pre- 
sent, however, are small and thus the variations in the 
coefficients would also be small. 

(9) A creep to failure phenomenon is observed that can be pre- 
dicted using a technique proposed by Crochet coupled with 
the modified Bingham model. 

(10) The proposed form of the loading-unloading equations for the 
modified Bingham model, Eqs. (3.3.3)', predicts the loading- 
unloading behavior of the adhesives reasonably well. 

(11) Initial results of relaxation tests at elevated temperatures 
indicate that the atdhesives follow the time-temperature 
superposition principle. 

(12) This study should give insight into rate and time dependent 
behavior of other adhesives and polymeric matrix materials 

- used in advanced composites as they are frequently similar 
in composition to the adhesives investigated herein. 

Future work may. be to compare the bulk properties of structural 
adhesives to those in the bonded state. An initial step would be to 
compare properties for single lap shear specimens to those for bulk 
shear specimens. Future endeavors in the area of loading-unloading 
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behavior should be very worthwhile. The effect of temperature on the 
rate properties should also prove to be an interesting area of future' 
research. Results from these tests could be correlated to those from 
the relaxation tests conducted at elevated temperatures herein. 
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